Kelvon-roton instability of vortex lines in dipolar Bose-Einstein condensates 
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The physics of vortex lines in dipolar condensates is studied. Due to the nonlocality of the 
dipolar interaction, the 3D character of the vortex plays a more important role in dipolar gases than 
in typical short-range interacting ones. In particular, the dipolar interaction significantly affects 
the stability of the transverse modes of the vortex line. Remarkably, in the presence of a periodic 
potential along the vortex line, a roton minimum may develop in the spectrum of transverse modes. 
We discuss the appropriate conditions at which this roton minimum may eventually lead to an 
instability of the straight vortex line, opening new scenarios for vortices in dipolar gases. 



Introduction.- When rotated at sufficiently large an- 
gular frequency, a superfluid develops vortex lines of zero 
density [H, [3] , around which the circulation is quantized 
due to the single-valued character of the corresponding 
wavefunction [3|. Quantized vortices constitute indeed 
one of the most important consequences of superfluidity, 
playing a fundamental role in various physical systems, 
as superconductors and superfluid Helium [E[ . Due to 
the enormous progress in the control and manipulation 
of ultra cold gases, Bose-Einstein condensates (BECs) of- 
fer an extraordinarily controllable system for the analysis 
of superfluidity, and in particular quantized vortices Q . 
Vortices and even vortex lattices have been created in 
BEC in a series of milestone experiments 0, B El • 

Vortex lines are indeed 3D structures, which, resem- 
bling strings, may present transverse excitations, which 
were studied for classical vortices by Lord Kelvin al- 
ready in 1880, and are therefore known as Kelvin modes 
[Toj | . These excitations have been also studied in the 
context of quantized vortices in superfluids by Pitaevskii 
ll| . Interestingly, the dispersion law for Kelvin modes at 
small wave vector k follows a characteristic dependence 
e(k) ~ — fc 2 ln/c£, where £ is the healing length. Kelvin 
modes play an important role in the physics of superfluid 
Helium [B|, [13] , and even of neutron stars [r| ■ Recently, 
Kelvin modes were experimentally observed in BEC [lj| . 

In spite of being extremely dilute, the properties of 
ultra cold g e.g. superfluidity, are fundamen- 

tally determined by the interatomic interactions. Up to 
very recently, typical experiments on ultra cold gases in- 
volved particles interacting dominantly via a short-range 
isotropic potential, which, due to the very low energies 
involved, is fully determined by the corresponding s-wave 
scattering length. However, recent experiments on cold 
molecules [Hf , Rydberg atoms [3] , and atoms with large 
magnetic moment [13], open a fascinating new research 
area, namely that of dipolar gases, for which the dipole- 
dipole interaction (DDI) plays a significant or even dom- 
inant role. The DDI is long-range and anisotropic (par- 
tially attractive), and leads to fundamentally new physics 
in condensates [3, E3, 0] , degenerated Fermi gases [21| , 
and strongly-correlated atomic systems [23]. It leads to 



the Einstein-de Haas effect in sphror condensates |23| . 
and may be employed for quantum computation [2~i| , and 
ultra cold chemistry [25j ] . Time-of-flight experiments in 
Chromium have allowed for the first observation ever of 
dipolar effects in quantum gases [2(5]. Remarkably, Fes- 
hbach resonances have been very recently employed to 
reduce the scattering length, and as a consequence the 
role of the dipolar interactions has been greatly enhanced 

Recently, the physics of rotating dipolar gases has at- 
tracted a growing interest. It has been shown that the 
critical angular frequency for vortex creation may be sig- 
nificantly affected by the DDI [2^]. In addition, dipolar 
gases under fast rotation develop vortex lattices, which 
due to the DDI may be severely distorted t 29], and even 
may change its configuration from the usual triangular 
Abrikosov lattice into other arrangements [3Cj . However, 
to the best of our knowledge, the previous analysis of vor- 
tices in dipolar BEC are constrained to situations where 
the 3D character of the vortex is unimportant. 

In this Letter, we analyze the role that the long-range 
character of the DDI plays in the physics of vortex lines 
in dipolar BEC. Due to this long-range character, differ- 
ent parts of the vortex line interact via DDI, and hence 
the 3D character of the vortices plays a much more im- 
portant role in dipolar gases than in usual short-range 
interacting ones. This has two main consequences. On 
one hand, for a fixed density the vortex core depends on 
the dimensionality of the problem, and can hence have in 
trapped gases a non-trivial dependence on the trap geom- 
etry. Additionally, depending on the dipole orientation, 
the DDI may significantly enhance or reduce the stiff- 
ness of the line against transverse excitations. Even more 
interestingly, under appropriate conditions, discussed in 
this Letter, the DDI may induce a minimum in the dis- 
persion law for the transverse modes, i.e. a Kelvon-roton 
spectrum. For sufficiently large DDI this minimum may 
reach zero energy, and hence the DDI may destabilize the 
straight-vortex configuration, opening the possibility for 
achieving other ground-state vortex configurations. 

Effective model.- In the following, we consider a dipo- 
lar BEC of particles with mass m and electric dipole d 
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(the results are equally valid for magnetic dipoles) ori- 
ented in the z-direction by a sufficiently large external 
field, and that hence interact via a dipole-dipole poten- 
tial: V d (f) = ad 2 (l — 3cos 2 #)/r 3 , where 9 is the an- 
gle formed by the vector joining the interacting particles 
and the dipole direction. The coefficient a can be tuned 
within the range — 1/2 < a < 1 by rotating the external 
field that orients the dipoles much faster than any other 
relevant time scale in the system [^H . At sufficiently low 
temperatures (and away from shape resonances (3^ |) the 
physics of the dipolar BEC is provided by a non-local 
non-linear Schrodinger equation (NLSE) of the form: 
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where g = Anh 2 aN /M , with a the s-wave scattering 
length (we consider a > 0). We assume the dipolar BEC 
in a ID optical lattice formed by counter-propagating 
lasers, V \(z) — sERsin 2 (qz), where Er — h 2 q 2 /2m is 
the recoil energy and q is the laser wave vector. 

In the tight-binding regime (sufficiently strong lat- 
tice) we can write ^(r, i) — J2j f( z — fyjOV'j ; (Pi *)> where 
b = ir/q, p — {x, y} and f(z) is the Wannier func- 
tion associated to the lowest energy band. Substituting 
this Ansatz in Eq. ([1]) we obtain a discrete non-linear 
Schrodinger equation. We may then return to a contin- 
uous equation, where the lattice is taken into account in 
an effective mass along the lattice direction and in the 



renormalization of the coupling constant [33L |3 
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where g — bg J f{z) i dz + g d C [11], with g d = aSird 2 /3, 
m* — h 2 /2b 2 J is the effective mass, and J = 
J f(z)[-(h 2 /2m)d 2 + V i(z)]f(z + b)dz. Note that the 
discreteness of the lattice is unimportant, i.e. Eq. ([2|) is 
valid, if k z <C 27r/6, being k z the z-momentum. In the 
following we use the convenient dimensionless parame- 
ter (3 = g d /g, that characterizes the strength of the DDI 
compared to the short-range interaction. It will be also 
useful in the following the Fourier transform of the DDI: 
V d (k) = g d [3cos 2 9 k - l]/2, with cos 2 9 k = k 2 /\k\ 2 . 

Homogeneous solution.- We consider first an homo- 
geneous solution of the form ^o(r,t) — y^exp[— i[it/h], 
where n denotes the condensate density, and fi — (g + 
V d (0))n is the chemical potential. From the correspond- 
ing Bogoliubov-de Gennes (BdG) equations one obtains 
that the energy e(k) corresponding to an excitation of 
wave number fulfills: e(k) 2 = i? kin (fc)[£? k i n (A;)+i?i n t(fc)] 
where E^i n (k) = hk 2 /2m+hk 2 /2m* is the kinetic energy, 




Figure 1: Dispersion e(q) as a function of q£, in the absence 
of an additional optical lattice, and various values of j3 = 
(solid), 1.0 (dashed), and -0.98 (dotted). 




Figure 2: Dispersion e(q) as a function of q£ for m/m* 
and p = -0.59. 



0.23, 



and Ei n t(k) = 2(g + V d (k))n. Stable phonons (i.e. exci- 
tations at low k) are only possible if Ei nt > for all k, 
i.e. if 2 + /?(3cos 2 6k - 1) > 0. If g d > phonons with 
k lying on the xy plane are unstable if (3 > 2, while for 
g d < phonons with k along z are unstable if (3 < — 1. 
Hence absolute phonon stability demands — 1 < (3 < 2. 

Vortex solution.- We consider at this point a conden- 
sate with a straight vortex line along the z-direction. The 
corresponding ground-state wavefunction is of the form 
^o(r,t) = (f>o(p) exjp(i(p) exp[— i/ji/H], where (p is the po- 
lar angle on the xy plane. The function </>o(p) fulfills 
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Mp) + 9\Mp)\ 2 Mp), (3) 



where D p = ^d p pd p and g = g~g d /2. Note that since 
is independent of z, the DDI just leads to a regularization 
of the contact interaction. The density of the vortex core 
is given by |<^ (p)| 2 , which goes to zero at p = 0, and be- 
comes equal to the bulk density n at distances larger than 
the corresponding healing length £ = h/^/mgn. Note 
that due to the regularization of the contact interaction, 
the size £ of the vortex core depends on the DDI. In par- 
ticular, this dependence is exactly the opposite as that 
expected for 2D vortices, since in 2D similar arguments 
provide g = g + g d . Hence, even for equal densities, and 
due to the long-range character of the DDI, the cores 
of a 2D vortex and a 3D vortex line can be remarkably 
different in a dipolar gas, differing significantly from the 
behavior of short-range interacting gases, where both 2D 
and 3D vortices would have the same core size. 
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The effects of the long-range character of the DDI be- 
come even more relevant in the physics of Kelvin modes. 
In order to analyze these modes, we consider a perturba- 
tion of the straight vortex solution of the form <3/(r, t) = 



*o (»*,*) + X{r, t) exp[i((p - fxt/K)], where [Til x(r,t) = 
Y.i\ u i{p) exp[i((pl + qz- et/h)] - vi{p)* exp[—i(<pl + qz - 
e*t/h)]]. Introducing this Ansatz into ((2]) and linearizing 
in x, one obtains the corresponding BdG equations: 
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with Fi(x,x') = Ii(x<)Ki(x>) , where and Ki 
are modified Bessel functions, and =max(x,x'), 
x < =min(x,a; / ). For every q we determine the lowest 
eigenenergy e, that provides the dispersion law discussed 
below. The first line at the rhs of Eqs. (j4|) and 
is exactly the same as that expected for a vortex in a 
short-range interacting BEC but with the regular- 
ized value g. The last term at the rhs of both equations is 
directly linked with the long-range character of the DDI 
and, as we show below, leads to novel phenomena in the 
physics of Kelvin modes (1 = 1) in dipolar BECs. 

In absence of DDI (or equivalently from Eqs. (j4} and 
([5]) without the last integral term) the dispersion law at 
low momenta (g£ -C 1) is provided by the well-known 
expression e(q) = —(h 2 q 2 /2m*)\n\[m/m*) 1 / 2 q(\. The 
integral terms of Eqs. (j4|) and ([5]) significantly modify the 
Kelvon spectrum in a different way depending whether 
P is positive or negative. In order to isolate the effect of 
the DDI on the core size with respect to the effect of the 
integral terms in Eqs. ([4]) and ([5]) we fix g and change the 
parameter fj which is proportional to the dipole-dipole 
coupling constant. Fig.[T]shows that for increasing (3 > 
the excitation energy clearly increases, i.e. the vortex 
line becomes stiffer against transverse modulations. 

In order to obtain an intuitive picture of why this is so, 
one may sketch the vortex core as a ID chain of dipolar 
holes. Dipolar holes interact in exactly the same way 
as dipolar particles, and hence maximally attract each 
other when aligned along the dipole direction, i.e. the 
^-direction. In this sense, the configuration of minimal 
dipolar energy is precisely that of a straight vortex line 
along z. A wiggling of the line produces a displacement 
of the dipolar holes to the side, and hence an increase 
of the dipolar energy. As a consequence, the DDI leads 
to an enhanced stiffness of the vortex line. From this 
intuitively transparent picture, we can easily understand 
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Figure 3: Stable/unstable regimes for straight vortex lines. 



that exactly the opposite occurs when j3 < 0. In that 
case, the dipolar holes maximally repel each other when 
aligned along the z-direction. Hence it is expected that 
the dipolar energy is minimized when departing from the 
straight vortex configuration. This results in a reduced 
energy of the excitations for j3 < ( Fig. [lj, i.e. it 
becomes easier to wiggle the vortex line. 

In principle, a sufficiently large DDI could even desta- 
bilize the straight vortex line. However, in the absence 
of an additional optical lattice along the z-direction, the 
destabilization would occur for values of /3 < — 1, i.e. in 
the regime of phonon instability in which the whole dipo- 
lar BEC is unstable against local collapses. An increase 
of the potential depth of the additional lattice leads to a 
reduction of the role of the kinetic energy term mq 2 jm* 
in Eqs. (g]) and ([5]) that enhances the effect of the dipo- 
lar interaction on the dispersion law. As a consequence, 
as shown in Fig. [21 in addition to the — k 2 ln(fc£) depen- 
dence at low momenta, a minimum eventually appears 
at intermediate momenta, i.e. the spectrum develops a 
Kelvon-roton character. Note that, typically, the disper- 
sion law presents a relatively abrupt change in its curva- 
ture for a given value of q (~ 1.4 in Fig. [2]). Interestingly 
this feature is given by an avoided level-crossing of the 
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lowest eigenvalues e(q). As a consequence the character 
of the lowest Bogoliubov modes changes at this value of 
q. For a sufficiently small (m/m*) cr the roton minimum 
eventually reaches zero energy, and the straight-vortex 
line becomes unstable against a novel type of instability. 
Fig.[3]shows as a function of (m/m*) and (i, the different 
regimes for the straight- vortex line. 

Roton minima occur in the dispersion law of superfluid 
Helium [lj , and have been also predicted in trapped dipo- 
lar gases under a strong ID and 2D confinement [13] • In 
those phonon-roton spectrum occurs related to 

excitations of the homogeneous BEC, and the eventual 
roton instability is linked to the appearance of density 
fluctuations (which are expected to lead to local collapses 
[35I]). Although the Kelvon-roton spectrum discussed in 
this Letter is also induced by the DDI, it fundamentally 
differs from the phonon-roton case. In this case, the roton 
minimum occurs due to the long-range interaction be- 
tween different positions in the vortex line. In addition, 
the eventual roton instability indicates that the straight 
vortex line ceases to be the ground-state configuration, 
opening the exciting possibility to achieve ground-state 
helicoidal vortex line configurations. The analysis of this 
possibility will be the subject of further investigations. 

Summary.- We have shown that, due to the nonlocal 
interaction between different parts of the vortex line, the 
3D character of the vortices is much more crucial in dipo- 
lar BECs than in short-range interacting ones. The size 
of the vortex core depends significantly on the system 
dimensionality, which indicates that, although we have 
restricted to homogeneous gases, vortex lines in trapped 
dipolar BECs may have a non-trivial dependence on the 
trap geometry. On the other hand, and depending on the 
dipole orientation, the DDI may significantly increase or 
decrease the stiffness of the line against Kelvin excita- 
tions. Interestingly, under appropriate conditions (ad- 
ditional optical lattice along the vortex line), a Kelvon- 
roton dispersion law occurs. For sufficiently large DDI 
a roton instability of the straight vortex line develops, 
which may preclude the appearance of a new type of he- 
licoidal ground-state configuration of the vortex line. 
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